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Abstract 



We formulate a real-space renormalization scheme that allows the study of the 
effects of bond randomness in the Heisenberg antiferromagnetic spin-1 chain. 
There are four types of bonds that appear during the renormalization flow. We 
implement numerically the decimation procedure. We give a detailed study 
of the probability distributions of all these bonds in the phases that occur 
when the strength of the disorder is varied. Approximate flow equations are 
obtained in the weak-disorder regime as well as in the strong disorder case 
where the physics is that of the random singlet phase. 
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I. INTRODUCTION 



The effect of quenclied impurities on tlie physics of one-dimensional spin systems is 
an important and unsolved problem. Many spin chains can be doped chemically and this 
creates some kind of disorder in the system. In addition the spin-1/2 chain is equivalent 
to a system of spinless fermions through the Jordan- Wigner transformation. This means 
that the problem of interacting spinless fermions in a disordered potential is equivalent to 
a random spin chain problem. There are not so many techniques that allow the study 
of these systems. The real-space renormalization group is prominent among them. Some 
time ago a pioneering study by Ma and Dasgupta |I| showed that the spin-1/2 Heisenberg 
antiferromagnetic chain with bond randomness is in a so-called random-singlet phase. In 
this phase the spins are locked into singlets that extend over arbitrarily long distances, in a 
pattern dictated by the bond distribution. It has been realized recently that the results of 
their renormalization procedure are in fact exact . This random- singlet phase may capture 
the physics of higher-dimensional disordered systems |^ . 

In the spin-1/2 case, the random-singlet phase appears for various kind of disorder and 
in a wide regions of the phase diagram when one adds XXZ anisotropy. This results from 
the study of the weak disorder regime by bosonizing the spin chain |Q. 

The spin-1 Heisenberg antiferromagnetic chain has a physics which is vastly different in 
the pure case. There is a gap to spin excitations and a finite spin-spin correlation length. 
These features can be best understood by consideration of a hidden topological order In 
fact the ground state of the spin-1 chain has a hidden long-range order that can be measured 
only by use of a nonlocal correlation function, the so-called string order parameter. It is 
a natural question to ask what happens to these peculiar features under the influence of 
disorder. In fact the original Ma-Dasgupta renormalization scheme requires a broad enough 
bond distribution to work 0. So more complex schemes have been proposed p|-p!0|. As a 
function of the disorder strength, it has been established that there is a phase transition 
between a low-disorder gapless phase with hidden order and a strong disorder phase which 
is the random-singlet phase of Ma and Dasgupta (gapless and no hidden order). 

In this paper, we give a detailed construction of a renormalization scheme suited to the 
study of the spin-1 chain. We generalize the Ma-Dasgupta decimation procedure by keeping 
more degrees of freedom. A brief account has been given in Letter form []TU| . Here we obtain 
explicit flow equations that are valid deep inside each of the phases that appear. We are 
able to follow the spin populations as a function of the renormalization scale as well as the 
evolution of distribution functions of the various kinds of bonds that appear. In section II, 
we define the renormalization scheme. In section HI, we study the weak-disorder phase of 
the spin-1 chain. Section IV contains our results for the strong disorder regime. The critical 
regime is studied in section V and section VI contains our conclusions. 



II. REAL SPACE RENORMALIZATION PROCEDURE FOR DISORDERED 
ANTIFERROMAGNETIC SPIN-1 CHAIN 

In this section, we explain how to obtain a real space renormalization scheme adequate 
to study the disordered antiferromagnetic spin-1 chain. 
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A. The Ma-Dasgupta real-space renormalization in the spin-1/2 case 



Ma and Dasgupta have introduced a real-space renormahzation procedure for the random 
antiferromagnetic spin-1/2 chain described by the Hamiltonian 



(2.1) 



where {Si} are quantum spin-1/2 operators and {Jj} are positive random variables dis- 
tributed with some probability distribution Po{J). Suppose that Ji is the largest coupling 
in the chain. The one-bond Hamiltonian, 



ho — J1S1.S2 — — 



' Si 



Ji 
2 



2 3 
~ 2 



admits two energy levels labeled by s = 0, 1 



Ji 
2 



Sis 



(2.2) 



(2.3) 



the level being (2s + 1) times degenerate : Cq = — represents the singlet, and ei = |Ji 
the triplet. At energies much lower than Ji, the spins 5*1 and S2 will therefore be frozen 
into the singlet state s = 0. The decimation procedure consists in eliminating the spins 5*1 
and 5*2, and in replacing the four spin segment Hamiltonian -^0,1.2, 3 involving the decimated 
spins 5*1 and S2 



-f^o,i,2,3 = ho + hi where hi = JqSq.Si + 725*2.5^ 
by the effective Hamiltonian for the remaining spins Sq and 5*3 



Hq!^/ — E'ns + JqSo-Ss 



(2.4) 



(2.5) 



which is meant to reproduce the four low-energy states of -f^o,i.2,3 which are separated from 
the other twelve states of -^0,1,2,3 by a big gap of order Ji. Using second order perturbation 
theory to treat hi gives 



E' 



0,3 



and 



Jo 



16J1 



J0J2 
'2J^ 



{Jl + Ji 



(2.6) 



(2.7) 



The same procedure may be iterated and successively applied to the new strongest bond of 
the chain. This defines a flow for the probability distribution of couplings P{J, ^l) where Q 
is the current strongest coupling [|I| 



dPjJ, Q) 



p{n, n) £ dJa £ dj, p{j,, n) p{j,, n)s(^j~ 



JgJb 

~2n 



(2.8) 



3 



This flow equation has to be supphed by some initial condition P(J, Qq). Fisher has shown 
that, for generic initial conditions, in the reduced variables F = In and 2; = ^ In (^j^ , 
the probability distribution R{z, F) of the variable z flows towards the unique fixed point 
R*iz) 

R{z, F) — y R*{z) = e{z) e"^ (2.9) 

r-+cxD 

where 6 is the Heaviside step function. This so-called Random Singlet Fixed Point corre- 
sponds to a power-law distribution in the original variables 

PV.f!) = 0(a - J) ^ (£)"'"'" where j-J^ (2.10) 

for which two typical bonds are typically much weaker than the strongest one Q. The ap- 
proximation involved in the use of perturbation theory to obtain the rule ( |2.7| ) therefore 
becomes better and better as the decimation proceeds, and the whole procedure is therefore 
completely consistent even if the initial distribution is not broad. The Ma-Dasgupta renor- 
malization scheme is moreover very appealing because it gives an interesting physical picture 
of the random spin-1/2 chain : at low energy, the chain is made of pairs of spins that are 
coupled together into singlets over arbitrarily long distances, the long singlets bonds being 
typically much weaker than the smaller ones. 



B. Renormalization of an AF bond between two spin-1 



The one-bond hamiltonian 



ho — J1S1.S2 — — 



SI 



Ji 
2 



(5^1 + S. 



^ \ 2 

2, 



(2.11) 



admits three energy levels labeled by s = 0, 1, 2 

Ji 



Y I'i' + 1) - 4] 



(2.12) 

2Ji represents the singlet, ei = — Ji 



the level being (2s + 1) times degenerate : cq = 
the triplet and 62 = Ji the quintuplet. 

In the Ma-Dasgupta procedure, there are only two levels, and "projecting onto the lowest 
level" is equivalent to "projecting out the highest level". Here these two possibilities are 
not equivalent. The first possibility has already been considered in refs.( [§], 0) where it is 
shown that the generalization of equation ( p.7| ) describing the effective coupling between Sq 
and 5*3 resulting from the projection onto the singlet formed by Si and S2, reads 



J'o 



4J0J2 
3 Ji 



(2.13) 



The coefficient | being bigger than 1, this rule is not automatically consistent : indeed, the 
inequalities Jo < Ji and J2 < Ji are not sufficient to imply that the new coupling Jq is 
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smaller than the decimated coupling Ji, in contrast with the rule ( |2.7| ) concerning spin- 1/2. 
This procedure can however be considered as qualitatively correct for very broad initial ran- 
domness, where the cases which would produce a new coupling Jq bigger than the decimated 
coupling Ji are statistically negligible. So the strongly disordered antiferromagnetic spin-1 
chains are described by the same random singlet fixed point already found in the study of 
disordered spin-1/2 chains. 

For weak initial randomness however, this naive procedure cannot be made consistent. 
We thus generalize the Ma-Dasgupta procedure with the interpretation of "projecting out 
the highest level" instead of "projecting onto the lowest level". More precisely for the 
antiferromagnetic bond described by the Hamiltonian h^, we project out the quintuplet 62 
but to keep the singlet cq and the triplet ei by replacing the two spin-1 Si and S2 by two 
spin-1/2 S[ and 5*2, and by replacing ho by the effective Hamiltonian 

K J 

= --f + JiS[.S', (2.14) 

The four spin segment Hamiltonian -f^o,i,2,3 containing the old spins Si and 5*2 

Ho,i,2,3 = ho + hi where hi = JoSo-Si + J2S2.S3 (2.15) 

has to be replaced by an effective Hamiltonian involving the spins So, S[, S'2 and S'3 

H&,s = hl" + hl^^. (2.16) 

If we use a first-order perturbation theory to treat hi, we find that the singlet of ho 
remains unchanged, whereas the degeneracy of the triplet is lifted by the perturbation hi. 
Using the Wigner-Eckart theorem for vectorial operators, we find more explicitly that the 
perturbation hi is equivalent to 

/if = QjoS^o + ^725^3) .(5^1 + 5^2) (2.17) 

We therefore have to choose the effective Hamiltonian 

hi" = JoSo.S'i + J2S'2.Ss (2.18) 

since it is equivalent at first-order perturbation theory, using again Wigner-Eckart theorem, 
to the hamiltonian 

[IJoSo + ^J2S3) {S'l + S^) (2.19) 

We have now enlarged the initial space since the chain now contains not only spin-1 but also 
spin-1/2. However it is possible to define a decimation procedure that is "closed" inside a 
particular set of spin chains as we will see in the following. 
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C. The real-space renormalization procedure 

We consider the enlarged set of spin chains described by the Hamiltonian 

= (2.20) 

i 

where the spin Si is a spin operator of size = | or Sj = 1, and where the couphngs {Jj} 

can be either positive or negative, but have to satisfy the following constraint : for any pair 
{i, j} such that i < j, the classical magnetization of the classical ground state of the segment 
must be smaller or equal to one in absolute value 



< 1 (2.21) 



where the quantity rrii j reads 



.J 

"n-l 



m. 



m-jp) 



p=i 



(2.22) 



n=i+l 

This condition for j = i + 1 gives immediately that there are exactly four types of bonds 

1) Link of type 1 : Ferromagnetic bond between two spin- 1/2 

2) Link of type 2 : Antiferromagnetic bond between two spin-1/2 

3) Link of type 3 : Antiferromagnetic bond between one spin-1 and one spin-1/2 

4) Link of type 4 : Antiferromagnetic bond between two spin-1 

Our decimation procedure is the following : 

To each bond (^Si, Si+i, Ji^ we associate the energy difference between the higher state 
and the lower state of the reduced Hamiltonian JiSi.Si^i 

Aj = — Jj if the bond i is of type 1 (2.23) 

Ai = Ji if the bond i is of type 2 (2.24) 

Aj = I Ji if the bond i is of type 3 (2.25) 

Aj = 3 Jj if the bond i is of type 4 (2.26) 

We pick up the bond {Si^^Si^^ corresponding to the strongest Aj of the chain. To 
define the renormalization rule for this bond, we again divide the four-spin Hamiltonian into 

Hio,h,i2,iz = ho + hi where ho = Ji^Si^Si^ and hi = Ji^Si^-Si^ + Ji^Si^.Si^ (2.27) 

and treat /ii as a perturbation of /iq to find the effective Hamiltonian replacing Hi^^i^^i^^i^ 
when the highest energy state of ho is removed. We have now to distinguish the four types 
of bonds 
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Rule 1) F bond between two spin-1/2 

The hamiltonian ho = Ji^Si^Si,^ admits two energy levels : the triplet ei = — and 

3|J- 

the singlet cq = The perturbation hi lifts the degeneracy of the triplet, and using 

Wigner-Eckart theorem, we find that hi is equivalent at first order of perturbation 
theory to 

hi = (^-JioSia + -</i2'Si3^ .{Si-^ + Si^) (2.28) 

To eliminate the singlet state and only keep the triplet state of ho, we remove the two 
spin-1/2 Si^ and St^ and replace them by a single spin-1 S'^_^, and we replace jg 

by 

<i.3 = + IJuA.A + IjJIAs (2.29) 



Rule 2) AF bond between two spin-1/2 

Here, we directly apply the Ma-Dasgupta procedure discussed in [II A| : we remove the 
two spin-1/2 Si-^ and Si^ and replace Hi^^i^^i^^i.^ by 



H. 



16 J, 



2Jn 



(2.30) 



Rule 3) AF bond between one spin-1 and one spin-1/2 

Suppose that = 1 and = |. The hamiltonian h^ = Ji^Si^Si^ admits two 

energy-levels : the doublet ei/2 = —Ji^ and the quadruplet 63/2 = "i^- At first order 
perturbation theory, Wigner-Eckart theorem gives that, within the subspace of the 
doublet s = 1, the perturbation hi is equivalent to 

hf = (aiJi^g + a2Ji^Si^ .{Sii + Si^) (2.31) 



where the constants ai and 0:2 read 

Si, + 1) - SiJs, 



ai 



1 + 



s(s + l) 



and 0:2 = 1 — ai = — 



1 



(2.32) 



The renormalization rule is therefore the following : we eliminate the spins S'j^ and 
5^2, and replace them by a single spin-1/2 5^^, and we replace ifjQ^ij_j2,i3 by the effective 
Hamiltonian 

eff 4 ^ 1 -> 

^io,i[,i3 ~ ~Jil + '^Jio^io^h ~ g ^^2 '^'ii 'S'jg . (2.33) 
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Rule 4) AF bond between two spin-1 

In this case we apply the rule explained at the beginning of this section (see eqs 
( p.l4| )-( pT8D ) : we replace the two spin-1 5*1 and ^2 by two spin-1/2 S*-^ and S[^, and 
we replace -f^o,i,2,3 by an effective Hamiltonian 



5Ji 



This renormalization procedure is entirely consistent from the point of view of the pro- 
gressive elimination of the highest energy degrees of freedom : it is easy to show that in the 
four cases of renormalization of a bond described above, all the energy scales Aj of the new 
bonds are always smaller than the energy scale Aji of the bond that we renormalize. 

It is also easy to check that this renormalization procedure is "closed" inside the set of 
spin chains defined by the condition ( |2.21| ) : if we apply this procedure to an initial chain 



belonging to this space, such as the random antiferromagnetic spin-1 chain we are interested 
in, the effective chain always belongs to this set of spin chains. In particular, spins higher 
than 1 cannot appear through this renormalization scheme. 

However, since this renormalization procedure is not purely based on complete decima- 
tion of bonds, it introduces correlations between bonds, so that it is impossible to write 
exact closed flow equations for the probability distributions of couplings, in contrast with 
the Ma-Dasgupta procedure. To study the properties of this renormalization scheme, we 
have therefore performed numerical simulations on spin-1 chains containing sites with pe- 
riodic boundary conditions (A^ = 2^^ for example), whose initial couplings Jj are distributed 
according to probability distributions of the following form 

P^( J) = i for 1 < J < 1 + d , and Pd{J) = elsewhere (2.35) 

The parameter d represents the strength of the initial disorder of the couplings. For a 
given number of sites N, and a given initial strength d of the disorder, we have numerically 
implemented the renormalization rules on a given number (typically 100) of initial indepen- 
dent samples, to compute averaged quantities over these different realizations of the initial 
disorder. It is convenient to use the variable : 

r = ln^, (2.36) 



where Q is the current strongest A (see eq. p.26|) and Qq the initial strongest A. We have 



studied the flow of the following quantities : the number A^(r) of effective spins = 1/2 
and S = 1 still present at scale F ; the proportion {A''(5'=i)(r)/iV(F)} of spins 5 = 1 among 
the effective spins at scale F ; the proportions Pi(T) = {Ni(T)/N(T)} of bonds of type 
i = 1, 2, 3, 4 at scale F; the probability distributions Pi{J, fl) of the coupling J at scale Q for 
the four types of bonds i = 1, 2, 3, 4. It is in fact more convenient to study the probability 
distributions Vi{x,T) of the reduced variable 

X = In (^) . (2.37) 
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where A (J) is defined as in ( p.26| ) 



A(J) = 
A(J) 
A(J) = 
A{J)- 



-J 
= J 

P 
3J 



for bonds of type 1 
for bonds of type 2 
for bonds of type 3 
for bonds of type 4 



(2.38) 
(2.39) 
(2.40) 
(2.41) 



so that the random variable x varies in (0, oo) for any type of bonds. 



III. THE WEAK DISORDER PHASE 

A. Numerical results 

In the weak disorder phase, we find that the number A^(r) of effective spins decays 
exponentially (see Fig |1|) 

N{T) oc e-"^'^)^ (3.1) 
r^oo 

where a{d) is a decreasing function of the disorder d that vanishes in the limit d —>■ d~ . As 
a consequence the magnetic susceptibility at temperature T can be computed by summing 
Curie laws for the free spins at scale Q = T. So we have : 

1 , , 

^^tT^- (3-2) 

The proportions Pi{T) of the four types of bonds reach a stationary regime characterized by 
(see Fig I) 

Pi(r)~o.25 p2(r)~o.75 P3(r)~o P4(r)~o (3.3) 

There are asymptotically only bonds of type 1 and bonds of type 2. This means in particular 
that there are only effective spin-1/2 in the chain, and no more spin-1. Since two bonds 
of type 1 cannot be neighbors according to the constraint ( |2.21| ), the even bonds and the 
odds bonds are not equivalent, as in the effective model of Hyman and Yang |^ : the 
"even" bonds are all antiferromagnetic, whereas the "odd" bonds are either ferromagnetic 
or antiferromagnetic with equal probability. 

It is necessary to introduce the probability distribution 'P|*'^"(x, F) for the couplings of 
the even bonds of type 2, and the probability distribution V2'^{x, F) for the couplings of the 
odd bonds of type 2. We find that 'P|''^"(x, F) becomes stationary for large enough F, and 
takes the form of an exponential distribution 

P|^""(x,F) e""^^ (3.4) 

where is independent of F, but depends on the value d of the disorder, and is numerically 
very close to the parameter a{d) characterizing the decay of A^(F) (|3.1| ). The probability 
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distributions Vi{x,T) and 7^2'^'^ {x,T) coincide (up to statistical fluctuations) and take the 
form of an exponential distribution (see Fig. |^) 



Pi(x, r) ~ '^(x, r) - a„(r) e-°°(^)" 

where the parameter ao(T) decays exponentially 

ao(r)oce-"'=^ 



(3.5) 



(3.6) 



As a consequence, for large enough F, the bond of the chain of highest A (corresponding 
to smallest x) that is chosen to be renormalized, is always an even bond of type 2. In the 
renormalization operation (2), this even bond disappear together with its two odd neighbors, 
and a new weak odd bond is produced. This explains why the distribution P|"^"'(x, F) for 
even bonds remains stationary, whereas the distribution of couplings of odd bonds becomes 
broader and broader in the variable x. This weak disorder phase is therefore the same as 
the "Haldane phase" found by Hyman and Yang in their effective model introduced in 0, 
and is very similar to the random dimer phase found in the study of random dimerized 
antiferromagnetic spin-1/2 chains |TT|] : in the asymptotic regime, the chain is made of a set 
of nearly uncoupled dimers. 



B. Approximate flow equations 

Assuming that the "even" bonds are all of type 2, that the "odd" bonds are either of 
type 1 or of type 2 with equal probability, and that the unique important process is the 
decimation of an even bond according to the rule (2) 



Sl 



S2 



S3 



1 « — 1 

2 S4 — 2 



Sl 



Ji 



J2 = n 



J3 



2n 



Si 



it is possible to write approximate flow equations for the probability distributions of the 
couplings are normalized according to 



1= dj p^'^'^un) = r dj p°'^\j,n) = r djPi{j,n) 

Jo Jo J-n 

It is convenient to introduce the normalized distribution of all odd bonds 

1 



p°'*'*( J, ^^) = 2 ( P2 U ^) + PiU^)] foi-n<j <n 



(3.7) 



(3.8) 



The approximate flow equations for the probability distributions P2'"'^"'{J, Q) and P°°'°'( J, Q) 
then read 



pf""(fi,fi) pf^"(j,fi) 



(3.9) 
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an 



-Pf""(fi,l])P°""^(J,l]) (3.10) 



+Pr"(^, ^) J^^ d.h Ji, d.h P°''\.h. n)6(^J- ^) (3.11) 

In the new variables F = In ^ and a; = In (Jj^ G [0, +oo), the flow equation for V2^^^{x) 
admits stationary solutions of exponential form 

P2"^<^«(a;) = 6-°-=^ (3.12) 

with undetermined constant Oe, in agreement with our numerical result ( p.4| ). With the last 
change of variables 

X — , z = ao{T)\n(^^ (3.13) 

the flow equation for the corresponding probability distributions Pi'^{z, T) and P2'^'^[z, T) 
admit the same stationary solution 

Pf\zJ) P^'^^zJ) — >e-' withao(r) oc e""^^ (3.14) 

r^oo r^oo 



where c^e is the number characterizing 1^2^^"^ {x) ( p.l2|) . We may also write the flow equation 



for the total number N(yt) of spins still present at scale Vt 

- — = -pr''{Si,n) iv(fi), (3.15) 

so that we obtain the following asymptotic behavior in the variable F : 

N{V) oc e-°^^. (3.16) 
r^oo 



IV. THE STRONG DISORDER PHASE 

A. Numerical results 

In the strong disorder phase d > dc, we find that the number A^(r) of effective spins 
decays as in the random singlet theory for the disordered antiferromagnetic spin-1/2 chain 
(see Fig |) : 

N{T)^ ^. (4.1) 
The magnetic susceptibility has thus the random singlet behaviour : 

Xoc (4.2) 

^ Tlog^T ^ ' 
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The proportions Pi{T) of the four types of bonds reach an asymptotic regime characterized 
by (see Fig |^) 

Pi(r) ~ P2(r) ~ e(r) p,{r) ~ 2e(r) p,{r) ~ i - 3e(r) (4.3) 

where e(r) slowly goes to as F ^ oo. This means that there is a sea of bonds of type 4, 
with sometimes defects of structure { bond of type 3, bond of type 2, bond of type 3 }. This 
defect structure is produced by the renormalization rule 4) for a bond of type 4 when its 
two neighbor bonds are also of type 4. The fact that there is no more bonds of type 1 in 
the asymptotic regime ( [4 .31 ) shows that defects are destroyed by the renormalization of the 
central bond of type 2 and not by the bonds of type 3 ; this means that for the probability 
distribution P4^{J,Q) at large enough fl, two typical couplings are much weaker than the 
bigger one. We indeed find that V4^{x, F) is an exponential distribution (see Fig ^ 

P4(a;,r) ~ a4(F)e-°^(^)^ (4.4) 

where the parameter l/a^lT) follows the random singlet behavior (see Fig. |^) 

^ F + Cst (4.5) 



a4(F) 



As a consequence, if a defect is produced at the renormalization energy scale Q, it survives 
until the energy scale j where it get decimated according to the rule (2), and the whole 
defect of structure {bond of type 3 , bond of type 2 , bond of type 3} entirely disappears 
to give one bond of type 4. Fig |^ shows indeed clearly that the probability distribution 
V2{x,T) tends to concentrate on the interval < x < ln3 as F increases. That has to be 
contrasted with the bonds of type 3, which are characterized by a distribution V^lx, F) that 
tends to coincide with Vi{x, F) for large enough F. 



B. Approximate flow equations phase 

Assuming that there is a sea of bonds of type 4, with sometimes defects of structure { 
bond of type 3, bond of type 2, bond of type 3 }, it is possible to write approximate flow 
equations for the probability distributions of the couplings normalized according to : 

1= / dj P2iJ,n)= [' dJP3iJ,n)= [' dj Pi{j,n). (4.6) 

Jo Jo Jo 

Assuming that the only two important renormalization processes are the production of the 
defect structure bond of type 3 , bond of type 2 , bond of type 3 by the renormalization rule 
4) for a bond of type 4 when its two neighbor bonds are also of type 4 



■So 



Sl 



S2 



S3 



Jo 



Jl - ^ 



J,. 



So 



S3 



Jo 



7 _ n 
Jl - 3 



J, 
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and the suppression of the defect structure by the decimation rule 2) 



So = 1 Si = I S2 = l S3 = 1 



So = 1 



S3 = l 



Jq Ji — J2 



Jo 



Jo J2 

2n 



we obtain the following approximate flow equations for the three probability distribu- 
tions : 



5{j-^]-P2{J,n) 



(4.7) 



on 3^3'; Nsin) 



p,{j,n) - p,{.j,n) 



(4.8) 



together with the flow equations for the number Ni{Q) of bonds of type « = 2, 3, 4 : 
dN2 



1^ = ^^^4 (f , ) N,{n) - P2{n, n) N2{n) 



(4.9) 
(4.10) 



(4.11) 



In 



P2{n,n) N2{n)-p,(^^,n^ N,{n) 



(4.12) 



so that the total number N{Q) = N2{^1) + N^i^) + N^IQ) of bonds evolves according to : 



dN 



-2P2in,n) N2in). 



(4.13) 



It is more convenient to write the flow equations for the probability distributions Vi{x, F) 
of the reduced variable a; = In 

(a?/))' ^^6^6 ^(J) is defined by (g^ so that the random 
variable x varies in (0, 00) for any type of bonds : 



^ - ^ ' +V2{0,T)V2ix,T) + -^V,iO,T) 



or 



dx 



N2{T) 



S{x -\n3) -V2{x,T) , 

(4.14) 



dT 



dV3(x,r) NAT) , ^ 
^ ' +2^f^ ^4(0,1) 



dx 



N2{T) 



V4{x-\n2,T)~V3{x,T) 



(4.15) 
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(9p4(x,r) dVi{x,T) 



dT 



dx 



+ 



iV2(r) 
"iV4(r) 



3^ 



(4.16) 



P2(o, r) dxi P3(xi, r) rfx2 ^3(^2, v) b[x-x^-x2- in -j . (4.17) 



Since the singular term containing the delta-function in ([4.14|) tends to develop a dis- 
continuity in ^2(3^, r) at a; = In 3, it is convenient to set : 



v^ix, r) = [1 - e(r)]^^i^^^^ + e(r)/2(x, r), 



(4.18) 



where < e(r) < 1 and /2(x, F) is a normalized probability distribution that is regular at 
X = ln3. Equation ( [4.14 ) will be satisfied if e(r) and /(a;, F) satisfy : 



e(r) = l-(ln3) ^^P4(0,r) 



de(r) 

a/(x,r) 9/(x,r) 



iV2(r) 
-e(r)[i-e(r)]/(o,r) 

+ /(o,r) /(x,r). 



9r dx 

Obvious stationary solutions for /(x, F) are simple exponentials : 

/(x,r)^ a/e-"/^ 

in which case e(r) vanishes exponentially : 

e(r) oc e""/^, 
r^oo 

so that V^ix^V) converges towards the stationary solution : 

^2(^) = l^^(ln3-x). 



This corresponds in the original variables to 

1 



for— < J <Q. 



(In 3) J 

We also obtain the following equation in the asymptotic regime 

1 



iV4(r)P4(o,r)~iV2(r) 



ln3' 



that we will use now to study the fiow equations for V^i^x, T) and V4,{x, T) 
With the last change of variables : 



(4.19) 

(4.20) 
(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 



(4.26) 



X — > z = a4(T)x, 



(4.27) 
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we find that the flow equation for the corresponding probability distributions P4,{z, T) and 
P3{z,r) admit the stationary solutions : 



where 



and Psiz, T) — > e 
r^oo 



1 



Q;4(r) oc - 
r^oo 1 



(4.28) 



(4.29) 



behavior of the total number A^(r) of spins (|4.13|) : 



as in the random singlet solution of Ma-Dasgupta. It is then easy to obtain the asymptotic 

(4.30) 



and the asymptotic behavior of the proportion e(r) (|4.3| ) of defects from ([4.26|) : 

N^iT) In 3 



(4.31) 



V. THE CRITICAL REGIME 

On Fig. ^ we have plotted the proportion ^'^^p'^''^'' of spins 5 = 1 among the effective 
spins for various values of the disorder this proportion flows towards in the weak disorder 
phase and to 1 in the strong disorder phase. Between these two attractive values, there is an 
unstable fixed point at dc — 5.75(5) where the proportion of spins S = 1 among the effective 
spins remains stationary at the intermediate value 0.315(5). The proportions Pi(r) of the 
four types of bonds reach a stationary state characterized by (see Fig ^) 

Pi(r)~o.i7, p2(r)~o.35, p3(r)~o.33, p4(r)~o.i5. (5.1) 

We find of course that the four probability distributions Vi{x,T) for i = 1,2,3,4 coincide 
up to statistical fluctuations (otherwise, the proportions Pi(r) would not remain stationary) 
and follow the exponential form (see Fig |TT|) : 

P,(x,r)~«e(r)e-"^(^)^ , (5.2) 

where the parameter 1/Q;c(r) (see Fig. |1^) follows the behavior of the effective model of 
Hyman and Yang : 

1 r 

~ - + Cst. (5.3) 



«c(r) 2 

The magnetic susceptibility is given by the effective number of free spins : 

Ycx ^^^T—. (5.4) 

^ Tlog^T ^ ^ 
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VI. CONCLUSION 



We have introduced a real-space renormalization scheme that allows the study of the 
spin-1 chain. Within this scheme we obtained a complete characterization of the weak- 
coupling phase, the critical regime and the strong-disorder phase. In all phases we were able 
to follow the spin populations and to obtain the probability distributions of the different 
types of bonds that appear under renormalization. It is only in the weak and strong coupling 
limit that we were able to obtain approximate analytical flow equations. 

The renormalization scheme that we used is an extension of the Ma-Dasgupta idea. These 
schemes have in common the fact that they are consistent for arbitrarily weak initial disorder. 
They do not create bonds stronger than the original decimated bond. In the spin-1/2 case, 
it is believed that this means that there is no critical disorder. In fact, this is suggested by 
bosonization: most bosonic forms of randomness give rise to relevant operators along the 
massless line of the pure system when the anisotropy is varied. The simplest assumption 0] 
is thus that the system flows immediately to the random-singlet phase (there is a region of 
stability of the spin liquid but this happens only for attractive enough interactions between 
the Jordan- Wigner fermions i.e. for negative enough anisotropy). 

However, this is not the case for the spin-1 chain. Here the Haldane gap is perturbatively 
insensitive to disorder as naively expected. This is known from bosonization studies of the 
spin-1/2 two-leg ladder as well as of the anisotropic spin-1 chain |[1^. So we may be in 



a situation with a first critical disorder strength corresponding to the vanishing of Haldane 
gap but which is unreachable by the real-space scheme. With increasing disorder there is 
then the second critical disorder strength for which the string order vanishes. This second 
transition is described by our renormalization scheme which is then asymptotically exact. 
Conversely, the bosonization methods are unable to follow the flow to strong coupling and 
thus are unable to describe even the weak-disorder phase captured by the real-space scheme. 
It may be also that there is nothing like a critical value of the disorder for the vanishing of 
the Haldane gap, if for example there are states of arbitrarily small energies in the gap as 
in the case of the Lifshitz tails in the localization problem. It remains to be seen if there is 
a single theoretical approach that is able to deal all known limiting cases. 
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FIG. 1. Linear-Log plot of the proportion of effective spins at scale F, for weak initial 
disorder d = 0.1, d = 0.5, d = 1 and d = 2 : this proportion decays exponentially (|3.1|). 
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FIG. 2. The proportions Pi(T) of the four types i = 1,2,3,4 of bonds at scale T, for weak 
initial disorder d = 0.1 : they reach the asymptotic regime ( |3.3[ ). 
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Probability distribution Ps^oddl^'H 

for weak initial disorder d=0.5 
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FIG. 3. Linear-Log plot of the probability distribution V2 (x, T) for F = 2, 3, 4, 5, 6, for weak 
initial disorder d = 0.5 : V2^'^{x, F) is well described by the exponential form ( ^.5| ) with a parameter 
ao(F) that is found to decay exponentially with F ( |3.6| ). 
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FIG. 4. Log-Log plot of the proportion of effective spins at scale F for strong initial 
disorder d = 100 : this proportion follows the power-law asymptotic behavior ([4.1|). 
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Proportions p^(r), p^{r), p^{r)/2, p^{r) 

for large initial disorder (d=100) 



^ Type 


i=4 


-o Type 


i=3 


^ Type 


i=2 


^ Type 


i=1 



10.0 



20.0 



30.0 

r 



40.0 



50.0 



60.0 



FIG. 5. The proportions Pi(T) of the four types i = 1,2,3,4 of bonds at scale F, for strong 
initial disorder d = 100 : they reach the asymptotic regime 
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Probability distribution P4(x,r) 



for large initial disorder (d=100) 
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FIG. 6. Linear-Log plot of the probability distribution Vi{x,T) for F = 4,8,12,16,20, for 
strong initial disorder d = 100 : 7^4(2;, F) is well described by the exponential form ( [4. 41) with a 
parameter 04 (F) plotted on Fig 0. 
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Probability distribution P2(x,r) 

for large initial disorder d=100 
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FIG. 8. Plot of the probability distribution 7^2(2:, F) for F = 8, 12, 16, 20, for a strong initial 
disorder d = 100 : V2{x,T) tends to concentrate on the interval < re < In 3 as explained in the 
text. 
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Proportion of spin-1 among the effective spins 

for initial disorder d=1 , 2, 4, 5.5, 6, 8, 1 6, 1 00 
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FIG. 9. Proportion of spins S = I among the effective spins at scale F for various values 
of the initial disorder d = 1,2,3,4,5.5,6,8, 16, 100 : this proportion flows towards in the weak 
disorder phase and to 1 in the strong disorder phase. Between these two attractive values, there 
is an unstable fixed point at dc — 5.75(5) where the proportion of spins S = 1 among the effective 
spins remains stationary at the intermediate value 0.315(5). 
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Proportions p^cr), p,(r), P3(r), p,(r) 

for critical initial disorder d =5.75 
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FIG. 10. The proportions Pi(T) of the four types i 
critical initial disorder dc = 5.75 (^). 



1,2,3,4 of bonds at scale T, for the 
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Probability distribution P2(x,r) 

for critical initial disorder d =5.75 

^ 

10 




0.0 10.0 20.0 30.0 40.0 50.0 60.0 

x=Ln(n/J) 



FIG. 11. Linear-Log plot of the probability distribution V2{x,T) for F 
critical initial disorder dc = 5.75 : 'P2{x,T) takes the exponential form 
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